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Magnetic resonance imaging with phase correction 
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The invention pertains to a magnetic resonance imaging system which 
involves correction of phase values of the magnetic resonance image. 



Magnetic resonance imaging with phase correction is known from the paper 'A 
new phase correction method in NMR imaging based on autocorrelation and histogram 
analysis ' by C.B. Ann and Z.H. Cho in IEEE Trans.MI 6(1987)32-36. 
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The known magnetic resonance imaging proposes to correct the phase of the 
complex valued pixels by first making an estimate of the first order phase error. This first 
order phase error is the contribution to the phase error in a complex valued pixel that is 
linearly dependent on the position of the pixel in the complex image. This first order phase 
error is calculated from the autocorrelation between adjacent pixels. On the basis of the 
estimate of the first order phase error a first-order phase error corrected image is derived. A 
histogram of phase values of the first-order phase error corrected image is formed From the 
location of one peak in this histogram or of peaks being separatedby 180° in this histogram, 
the zero-order phase error is estimated. This zero-order phase error is the overall contribution 
that is independent of the position of the position of the pixel in the complex image. 



An object of the invention is to provided an magnetic resonance imaging 
method in which phase errors having a complicated nature are corrected more accurately. 

This object is achieved by a magnetic resonance imaging method according to 
25 the invention comprising a reconstruction unit arranged to 

- reconstruct a complex image of complex valued pixels from magnetic resonance 

signals 

- compute a distribution of phase values of the complex image 
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- apply a phase correction to the complex image to form a corrected magnetic 
resonance image and 

- control the phase correction on the basis of the distribution of phase values of the 
complex image. 

According to the invention the phase correction is controlled on the basis of 
the distribution of me phase values of the complex image. It appears that the distribution of 
phase values is quite sensitive for the way the phase error depends on the position in the 
complex valued image. Consequently, the distribution of phase values is an accurate indicator 
for the kind of phase corrections to be effective and also for the accuracy of the result 
achieved by the phase correction actually applied. These and other aspects of the invention 
will be further elaborated with reference to the embodiments defined in the dependent 
Claims. 

Preferably, the distribution of phase values in the complex image is 
represented by a histogram of phase values. Generation of the histogram from the complex 
values of the pixels of the complex image requires very low computational effort and takes 
only a small memory capacity. Further, the features of the histogram appear to be very 
sensitive to the way the phase error varies over the complex valued image. Moreover, the 
features of the histogram also accurately indicate the accuracy achieved in correction of the 
phase error. 

The use of test functions applied to the histogram enhances the sensitivity for 
particular patterns of variation of the phase error over the complex valued image. Further, the 
test function provides an easy indication of the accuracy of the correction in terms of a simple 
real valued number. Particularly suitable test functions are chosen from the class of functions 
/that satisfy: f(h)>h (h>l). These test functions enhance peaks in the histograms 
relative to broader structures. Accurate correction of the phase errors leads to peaks in the 
histogram and removes broad structures. As the test functions of this aspect of the invention 
favour peaks in the histogram, these test functions provide an accurate indication of the 
accuracy of the correction of the phase errors. Particularly good results are obtained when the 
histogram power function is employed as the test function. The histogram power function is 
defined as P(H) = £hj . Here JJdenotes the histogram and hj is the number of phase- 
values in the j-th bin of the histogram H. The histogram power function strongly emphasises 
histograms having a few strong peaks and returns low values for histograms having mainly 
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broad features. Thus the histogram power function provides an effective test which indicates 
if the phase correction is accurate. 

In a further implementation of the invention polynomial phase corrections are 
carried-out Such polynomial phase corrections require low computational effort and require 
only a short computation time. Accordingly, several polynomial phase corrections can be 
tried within a short time and for the respective phase corrections the histogram is constructed 
and the histogram power function is computed. Then, the most effective correction is selected 
on the basis of the minimum histogram power function value. The respective polynomial 
phase corrections are easily effected just by adapting the polynomial coefficients of the 
polynomial phase corrections and re-iterating the phase correction with the adjusted 
polynomial coefficients. Preferably, the polynomial coefficients are adjusted so as to 
optimise me value of the test function of the histogram. These polynomial coefficients for 
which the test function of fee histogram acquires its optimum value correspond to the most 
effective polynomial phase correction. 
15 The computational efficiency of the phase correction is further improved in 

that the phase correction is adjusted by resetting the polynomial coefficients on the basis of a 
trial and improve algorithm so arrive at a set of polynomial coefficients that correspond to a 
phase distribution of which the test function of the histogram has it optimum value, notably 
the histogram power function, is minimum . 
20 The corrected magnetic resonance image is practically void of artefacts caused 

by phase errors in the magnetic resonance signals. Accordingly, the corrected magnetic 
resonance image has a high diagnostic value in that the anatomy of the patient to be 
examined is imaged faithfully, even for small details with low contrast or low contrast- 
resolution/ 

25 The invention also relates to a magnetic resonance imaging method as defined 

in Claim 10. The magnetic resonance imaging method of the invention achieves accurate 
correction of phase errors in the complex valued image. The invention further relates to a 
computer programme as defined in Claim 1 1 . The computer programme of the invention can 
be loaded into the reconstruction unit of the magnetic resonance imaging system. 
30 Accordingly the reconstruction unit is able to perform accurate correction of phase errors in 
the complex valued image feat is reconstructed from the magnetic resonance signals 
generated by the magnetic resonance imaging system. 
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These and other aspects of the invention will be elucidated with reference to 
the embodiments described hereinafter and with reference to the accompanying drawing 



wherein 



Figure 1 shows a flow chart of an implementation of the method of the 

invention and 

Figure 2 shows a schematic representation of a magnetic resonance imaging 
system in which the invention is employed. 



Figure 1 shows a flow chart of an implementation of the method of the 
invention. The magnetic resonance imaging method of the invention involves acquisition 
(mr-acq)of magnetic resonance signals. From the acquired magnetic resonance signals, the 
complex image is reconstructed (tm-recon). For individual complex valued pixels of the 
complex image, the phase is computed. Then the phase-values for the pixels are binned into 
th histogram. That is, phase-values in pre-determined phase-intervals, the bins, are counted 
and the numbers of phase values in respective bins form the histogram H. From the histogram 
the histogram power function P(H) = £/jf is computed and on the basis of the valueof 



the histogram power function a polynomial phase correction is selected and applied to the 
complex image. In an iterative way, again the histogram and its histogram power function are 
computed and a new polynomial phase correction is selected and applied to the current 
complex image. The iteration is ended when the histogram power function has its minimum 
value which indicates that the histogram contains a few strong peaks which indicates that the 
current phase correction is adequate. Optionally, a final phase correction can be performed on 
the basis of the histogram having the minimum histogram power function value. For 
example, the optional final phase correction may involve an accurate water-fat shift 
correction. This water-fat shift correction is based on the separation of two main peaks in the 
histogram. 

Figure 2 shows a schematic representation of a magnetic resonance imaging 
system in which the invention is employed. The magnetic resonance imaging system includes 
a set of main coils 10 whereby the steady, uniform magnetic field is generated. The main 
coils are constructed, for example in such a manner that they enclose a tunnel-shaped 
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examination space. The patient to be examined is slid into this tunnel-shaped examination 
space. The magnetic resonance imaging system also includes a number of gradient coils 1 1, 
12 whereby magnetic fields exhibiting spatial variations, notably in the form of temporary 
gradients in individual directions, are generated so as to be superposed on the uniform 
5 magnetic field. The gradient coils 11, 12 are connected to a controllable power supply unit 
21. The gradient coils 11, 12 are energised by application of an electric current by means of 
the power supply unit 21 . The strength, direction and duration of the gradients are controlled 
by control of the power supply unit The magnetic resonance imaging system also includes 
transmission and receiving coils 13, 16 for generating the RF excitation pulses and for 
10 picking up the magnetic resonance signals, respectively. The transmission coil 13 is 

preferably constructed as abody coil 13 whereby (apart of) the object to be examined can be 
enclosed. The body coil is usually arranged in the magnetic resonance imaging system in 
such a manner that the patient 30 to be examined is enclosed by the body coil 13 when he or 
she is arranged in the magnetic resonance imaging system. The body coil 13 acts as a 
15 transmission antenna for the transmission of the RF excitation pulses and RF refocusing 

pulses. Preferably, the body coil 13 involves a spatially uniform intensity distribution of the 
transmitted RF pulses (RFS). The same coil or antenna is usually used alternately as the 
transmission coil and the receiving coil. Furthermore, me transmission and receiving coil is 
usually shaped as a coil, but other geometries where the transmission and receiving coil acts 
20 as a transmission and receiving antenna for RF electromagnetic signals are also feasible. The 
transmission and receiving coil 13 is connected to an electronic transmission and receiving 
circuit 15. 

It is to be noted mat it is alternatively possible to use separate receiving coils 
16. Forexample,surfececoilsl6canbeusedasreceiving coils. Such surface coils have a 
25 high sensitivity in a comparatively small volume. The transmission coils, such as foe surface 
coils, are connected to a demodulator 24 and foe received magnetic resonance signals (MS) 
are demodulated by means of foe demodulator 24. The demodulated magnetic resonance 
signals (DMS) are applied to a reconstruction unit The receiving coil is connected to a 
preamplifier 23. The preamplifier 23 amplifies foeRF resonance signal (MS) received by foe 
30 receiving coil 16 and foe amplified RF resonance signal is applied to a demodulator 24. The 
demodulator 24 demodulates foe amplified RF resonance signal. The demodulated resonance 
signal contains foe actual information concerning foe local spin densities in foe part of foe 
object to be imaged Furfoermore, foe transmission and receiving circuit 1 5 is connected to a 
modulator 22. The modulator 22 and foe transmission and receiving circuit 1 5 activate foe 



r±llNUHU372KP.P 



6 14.04.2003 
transmission coil 13 so as to transmit the RF excitation and refocusing pulses. The 
reconstruction unit derives one or more image signals from the demodulated magnetic 
resonance signals (DMS), which image signals represent the image information of the 
imaged part of the object to be examined The reconstruction unit 25 in practice is 
constructed preferably as a digital image processing unit 25 which is programmed so as to 
derive from the demodulated magnetic resonance signals the image signals which represent 
the image information of the part of the object to be imaged. The signal on the output of the 
reconstruction monitor 26, so that the monitor can display the magnetic resonance image. It is 
alternatively possible to store the signal from the reconstruction unit 25 in a buffer unit 27 
while awaiting further processing. 

The reconstruction unit of the magnetic resonance imaging system according 
to the invention is also provided with a processing unit, for example in the form of a 
computer which includes a (microprocessor. The processing unit of the reconstruction unit 
is arranged to carry out the phase correction to the complex valued image on me basis of the 
distribution of phase values in the complex valued image, e.g. as represented by the 
histogram To this end, the computer program according to the invention is loaded, for 
example, into the processing unit of the reconstruction unit 25. 
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CLAIMS: 



10 



1. A magnetic resonance imaging system comprising a 

reconstruction unit arranged to 

- reconstruct a complex image of complex valued pixels from magnetic resonance 

signals 

_ compute a distribution of phase values of the complex image 

- apply a phase correction to the complex image to form a corrected magnetic 
resonance image and 

- control the phase correction on the basis of the distribution of phase values of the 
complex image. 

2. A magnetic resonance imaging system as claimed in Claim 1, wherein Ihe 

distribution of phase values of 1he complex image is represented by a histogram of the phase 
values of Ihe complex image. 

15 3- A magnetic resonance imaging method as claimed in Claim 2, wherein the 

phase correction is controlled on the basis of a test function of the histogram. 

4. A magnetic resonance imaging system as claimed in Claim 3 wherein the test 
function of the histogram discriminates, in particular enhances peaks, peaks in the histogram 

20 from broader distributions. 

5. A magnetic resonance imaging system as claimed in Claim 4, wherein the test 
function is formed by the histogram power function. 

25 6. A magnetic resonance imaging system as claimed in Claim 1, wherein the 

reconstruction unit is arranged to make the phase correction on the basis of apolynomial 
phase correction, said polynomial being represented by its polynomial coefficients. 
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7. A magnetic resonance imaging method as claimed in Claim 6, wherein the 
reconstruction unit is arranged to control the phase correction by adjusting polynomial 
coefficients of the polynomial phase correction. 

8. A magnetic resonance imaging system as claimed in Claim 3, wherein the 
reconstruction unit is arranged 

- to make the phase correction on the basis of a polynomial phase correction and 

- control the phase correction by adjusting polynomial coefficients of the polynomial 
phase correction so as to optimise the test function. 

9. A magnetic resonance imaging system as claimed in Claim 8, wherein the 
polynomial coefficients are adjusted by way of a trial and improve algorithm controlled On 
the basis of the test junction. 



A magnetic resonance imaging method wherein 
a complex image of complex valued pixels is reconstructed from magnetic 
resonance signals 

a distribution of phase values of the complex image is computed 

a phase correction is applied to the complex image to form a corrected magnetic 

resonance image and 

the phase correction is controlled on the basis of the distribution of phase values of 
the complex image. 

A computer programme comprising instructions to 
compute a distribution of phase values of a complex image 
apply a phase correction to the complex image to form a corrected magnetic 
resonance image and 

control the phase correction on the basis of the distribution of phase values of the 
complex image. 
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ABSTRACT: 



A magnetic resonance imaging system comprises reconstruction unit that is 
arranged to reconstruct a complex image of complex valued pixels from magnetic resonance 
signals. A compute a distribution of phase values of the complex image and to apply a phase 
correction to the complex image to form a corrected magnetic resonance image. The phase 
correction is controlled on the basis of the distribution of phase values of the complex image. 
Notably, toe histogram power function is an effective indicator of the accuracy of the phase 
correction. 
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